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Improved State Estimation in Quadrotor MAVs: A 
Novel Drift-Free Velocity Estimator 

Dinuka Abeywardena, Sarath Kodagoda, Gamini Dissanayake, and Rohan Munasinghe 


Abstract —This paper describes the synthesis and evaluation 
of a novel state estimator for a Quadrotor Micro Aerial Vehicle. 
Dynamic equations which relate acceleration, attitude and the 
aero-dynamic propeller drag are encapsulated in an extended 
Kalman filter framework for estimating the velocity and the 
attitude of the quadrotor. It is demonstrated that exploiting 
the relationship between the body frame accelerations and 
velocities, due to blade flapping, enables drift free estimation of 
lateral and longitudinal components of body frame translational 
velocity along with improvements to roll and pitch components 
of body attitude estimations. Real world data sets gathered using 
a commercial off-the-shelf quadrotor platform, together with 
ground truth data from a Vlcon system, are used to evaluate 
the effectiveness of the proposed algorithm. 

Introduction 

Quadrotor Micro Aerial Vehicles (MAV) are simple robotic 
platforms to construct. In its’ basic form, it is no more than 
two counter rotating propeller pairs attached symmetrically 
to a rigid cross-like frame, along with the means to control 
the speed of each individual propeller. This symmetric design 
is what has enabled the quadrotor to become a simple yet 
powerful vertical take-off and landing aerial platform that is 
popular among the robotics community. 

With this simplicity comes the burden of controlling motion 
in 3D space with the use of just four actuators. Under¬ 
actuated and coupled dynamics of the quadrotor make it nearly 
impossible for a human pilot to gain control of it, unless a well 
tuned control system is in place. Such a control system is also 
vital if autonomy is a goal, as is the case with most MAVs. 
Estimates of controlled states and their derivatives are essential 
for any control system, and where those estimates are accurate 
and frequent in time, it has been demonstrated that quadrotors 
have extreme maneuverability and agility JT]. 

However, MAVs are - by design - limited in their payload 
capacity and with those limitation, obtaining accurate and fast 
state estimates becomes a challenge. For example, MEMS 
inertial sensors can provide fast but coarse state estimates 
0, while exteroceptive sensors such as lasers and cameras 131 
render more accurate state estimates albeit at a slower rate. 
Attempts to merge these two sensing domains are frequent in 
MAV literature H, 0 and an application of similar ideas to 
quadrotors was presented in ||6]. 

One aspect common to most MAV state estimators is their 
use of inertial sensors. Typically gyroscopes, accelerometers 
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and magnetometers are used for the purpose of attitude es¬ 
timation Q. Based on a long history of research in inertial 
navigation systems, sensor fusion algorithms usually employed 
for this task make use of the equations of motion of the 
sensing unit in three-dimensional space. The main advantage 
of this approach is that these generic estimators are specific 
only to the sensor package geometry and as such can be 
used independently of the platform on which the sensors are 
mounted. However, they fail to exploit the dynamics of the 
vehicle under consideration in the estimation process, leading 
to a potentially sub-optimal result. The value of using specific 
dynamic characteristics of the vehicle has been reported in the 
case of land vehicles 10 and air vehicles®. Similarly, in this 
paper we demonstrate that the influence of blade flapping in 
a quadrotor leads to a set of dynamic equations that can aid 
state estimation using inertial sensors. 

The rest of the paper is arranged as follows. First, some 
background on quadrotor state estimators and a discussion 
on what motivated us to look beyond the state-of-the-art is 
presented. We will then briefly present the quadrotor dy¬ 
namic equations that are of interest to the state estimation 
process. After highlighting the shortcomings of the generic 
design, a novel state estimator design is presented along 
with experimental results which demonstrate the accuracy and 
consistency of estimates. The article is conclude by exploring 
the implications of the novel algorithm. 

Background and Motivation 

MAV attitude estimators that fuse gyroscope and accelerom¬ 
eter measurement using generic algorithms are frequently re¬ 
ported in literature 1101,0,0. In a nutshell, these algorithms 
operate by fusing measurements of a triad of body mounted 
gyroscopes and accelerometers. Gyroscope measurements are 
a source of high frequency attitude rate information, but they 
alone are not sufficient for drift free attitude estimation due 
to bias and various other forms of noise present in a typical 
low cost sensor. Attitude estimators for MAVs overcome this 
issue by assuming that accelerometers predominantly measure 
gravitational acceleration and are thus capable of providing 
low frequency information about MAV orientation with re¬ 
spect to gravity. Clearly, when the vehicle accelerations are 
significant, as in the case of quadrotor, this assumption does 
not hold im. Furthermore, such estimators are incapable of 
drift free velocity estimation, as they can only be generated by 
integrating noisy accelerometer measurements. To complicate 
the matters even further, accelerometer measurements need 
to be compensated for gravity before this integration, and 
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such a compensation requires an accurate attitude estimate. 
As mentioned before, one promising way to overcome these 
deficiencies is to examine the behaviour of the MAV in 
question, in-order to identify suitable characteristics that would 
assist the estimation process. 

Martin et. al. ifT^ have analysed the behaviour of a quadro- 
tor MAV in detail and also presented equations describ¬ 
ing measurements of an accelerometer mounted on-board a 
quadrotor. Their results motivated us to reformulate the state 
estimators for quadrotors and to redesign them considering the 
true sensor behaviour as opposed to conventional vehicle inde¬ 
pendent assumptions. In addition to improving the accuracy of 
the attitude estimate, the design presented here provides a drift 
free estimate of the horizontal components of translational ve¬ 
locity of the quadrotor. Recently, a similar idea was presented 
in ||T 3 ]| where two separate non-linear complementary filters 
were utilised to estimate attitude and velocity of a quadrotor 
MAV. The filter formulation presented in this paper is different 
from Ha and we also present experimental results validating 
the concept. The velocity estimates thus derived are of critical 
importance to control and navigational tasks of a quadrotor, 
as will be discussed in concluding remarks. 

Quadrotors: What Makes Them Unique? 

A thorough derivation and analysis of the quadrotor dy¬ 
namics can be found in HI and Ha. Rather than re¬ 
iterating the derivation, here we aim to briefly summarise the 
important equations and to provide an intuitive description of 
the most salient features of the dynamic behaviour that makes 
quadrotors a unique MAV. 

Let {E} be the earth fixed inertial frame, and a vector 
[ X y z Y" denote the position of the centre of mass of 
the quadrotor expressed in {E} (See Fig. [^. Let {B} = 
[ bi 62 ^3 be a body fixed frame positioned at the centre 

of mass of the quadrotor. 



Fig. 1. Coordinate frame definitions for the quadrotor dynamic model 

The orientation of { 5 } with respect to {E} is defined using 
a cumulative rotation of Euler angles tjj (Yaw) , 6 (Pitch) and 
(p (Roll) in that order, around 63, 62 and bi, respectively. R 
is defined as the rotational transformation matrix from {B} 
to {E}. The kinematic equation relating the instantaneous 


angular velocity U = [ ujy \ of {B} with respect to 
{E}, to Euler rates can be expressed as: 
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The equation describing the evolution of translational mo¬ 
tion of the quadrotor as derived in ca is of special interest 
to the estimator design that will be presented in following 
sections. 

4 4 

mV = mg — kr ^ — Ai ^ w^V ( 2 ) 

i=l i=l 

where 

V = Velocity of {B} as observed from an inertial frame 
g = gravity vector 

kr = thrust coefficient of propellers 

Al = a positive coefficient known as rotor drag coefficient 
uJi = rotational velocity of rotor, i G { 1 , 2 , 3 , 4 } 

V = projection of V on to the propeller plane 
m = mass of the quadrotor 

Equation Q sheds light on two key aspects of the quadrotor. 
Eirst and the most obvious is the fact that the thrust force is 
perpendicular to the propeller plane, and thus has no effect 
on motion along that plane. Secondly and more importantly, 
we see the presence of a force which is proportional to 
the translational velocity of the quadrotor. For an intuitive 
description of this force, we refer readers to Fig. which 
shows a cross section of a quadrotor in flight, and provide 
below a simplified explanation of the origin of this force. 



Fig. 2. Schematic of a quadrotor immediately after tilting sideways, but 
before it starts moving. ET is the summation of propeller thrusts and 
corresponds to the second term in 

Fig. shows a quadrotor in a hypothetical state where it has 
tilted sideways to initiate a translation in a horizontal direction 
but immediately before it gains any translational motion. At 
this point thrust from propellers and gravity are the only forces 
acting on our simplified quadrotor model. As is obvious, in 
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this particular state, thrust force generated from the propellers 
is perpendicular to the propeller plane. 



Fig. 3. Quadrotor, after tilting, starts moving sideways. fi and /2 are the 
orthogonal components of ET 


direction of body 
motion 



position of max 
flapping angle 


the amount of blade flapping is dependent on the translational 
velocity of the quadrotor, the component of thrust force along 
the body plane is also a function of that velocity. The last 
term in Q models the impact of this component of the thrust 
force on the translational motion of the quadrotor. If one is 
to place an accelerometer on-board the quadrotor, with its’ 
sensing axis parallel to the propeller plane, that accelerometer 
will measure a force that is roughly proportional to the velocity 
of the quadrotor along the same axis. In fact, in the next 
section, it is shown that this is the only significant force that 
the said accelerometer will sense. (Interestingly, Q ignores the 
aerodynamic drag experienced by a body moving through air, 
which is usually a function of the square of the velocity. This 
can be justified for quadrotor MAVs that move at relatively low 
speeds.) This is the unique characteristic of quadrotor MAVs 
that will later be exploited to the benefit of the state estimator. 

To conclude this section, we re-write (j^ using ^V (i.e. 
V in {B} frame) to facilitate the estimator design. After 
neglecting the second order terms that appear due to co¬ 
ordinate frame transformation, the first two components of 
^V S Vy,^ Vz} can be written as 


b ■ ■ a b 

Vx ~ —9 sin t) - Vx 

m 

b ' /i ‘ A b 

Vy ~ g cos Cl sin 0- v. 

m 


where 


(3) 
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fci = Ai ^ uj, 

i=l 


In what follows, we assume that ki is a positive constant 
considering the fact that the summation of propeller rotational 
rates are fairly constant during smooth flight. 


position of max 
flapping velocity 


Fig. 4. As the propeller blades rotate, flapping is determined by their position 
with respect to the direction of motion of the propeller as a whole. 

As stated, the state depicted in Fig. is hypothetical in 
the sense that even the slightest tilt of the quadrotor will 
induce translational motion. Fig. shows a more realistic 
situation in which quadrotor now moves right with a non¬ 
zero velocity. For a propeller with two blades, we can now 
identify a retreating and an advancing blade, as shown in 
blue and green respectively in Fig. The velocity of the 
advancing blade with respect to free air is higher than that 
of the retreating blade, due to the translational velocity of the 
whole quadrotor. This creates a force imbalance between the 
two blades of the same propeller and thus causes the blades 
to flap up and down as they rotate. Blade flapping forces the 
propeller to rotate out of plane and the flapping angle of a 
blade is at a maximum just before it transitions from advancing 
state to retreating state or vice versa. As shown in Fig. blade 
flapping causes the thrust force of the propeller to be tilted in 
a direction which opposes the motion of the quadrotor. As 


Inertial Sensors in Quadrotors 

This article is concerned with the quadrotor state estimators 
based on inertial sensors and specifically with accelerometers 
and gyroscopes. For simplicity, we assume that a triad of 
accelerometers and gyroscopes are mounted at the centre of 
mass of the quadrotor body. For both types of sensors we 
adhere to standard MEMS error models HD. 

Gyroscopes measure the instantaneous rotational rate of the 
body with respect to the inertial frame, and their measurements 
can be modelled independently of the equations of motion of 
the moving platform to which they are attached. 

gi = Pgi + Wgi (4) 

i^gi — i^gi ^/3gi ( 5 ) 

Tgi 

where jSgi is the bias of gyroscope and Tgi is the time 
constant of gyroscope bias. Wgt and wpgi are zero mean 
White Gaussian Noise (WGN) terms. 

In contrast, accelerometers measure a combination of in¬ 
ertial and gravitational acceleration, and their measurements 
can be expressed using the equations of motion governing 
the body they are mounted on. Perhaps one of the best 
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example of the value of this strategy is the case of a triad 
of accelerometers mounted on a quadrotor platform. Denoting 
by di the acceleration that would be measured by an ideal 
accelerometer, we combine the accelerometer measurement 
model with Q to arrive at: 

a = V-g 

4 4 

= ( 6 ) 

i=l i=l 

Equation (|^ describing the readings obtained from an on¬ 
board triad of accelerometer is unique to quadrotors and is 
of critical importance to a state estimator in that context. As 
stated in the previous section, equation shows that the 
accelerometers along bi and b 2 coordinate axes are only 
sensitive to a force which is dependant on the projection 
of the quadrotor translational velocity on to bi, b 2 plane. 
Furthermore, the component of the gravitational acceleration 
in the body frame (which is typically large compared to inertial 
accelerations of slow moving vehicles) no longer influences 
the accelerometer measurement. In the following section, we 
will exploit this unique property to design a better state 
estimator for quadrotors. 

Estimator Design 

The goal here is to design a state estimator for the quadrotor, 
giving due regard to the dynamic and kinematic equations 
presented in the previous sections. For this, we propose a six 
state. Extended Kalman Filter (EKE) based state estimator. 
The filter states are: 


(f) — Roll angle in current orientation estimate 
9 — Pitch angle in current orientation estimate 
Pgx — Bias in X axis gyroscope 
Pgy — Bias in Y axis gyroscope 

^Vx — X velocity component of quadrotor in body frame 
^Vy — Y velocity component of quadrotor in body frame 

Process Model 

Equations ([T]i, ([^ - Q form the EKE process equation. 
Out of the three Euler angles we can only estimate (j) and 9 
as the equations are expressed in a form independent of the 
yaw angle ijj. 

4'= i9x- Pgx + Wgx) + tan 9 cos - /Sgz) I 

+ tan 9 sin (j){gy — jdgy + Wgy) \ (!) 

9 = cos 4>{gy - 4igy + W gy) - SlO ) ) 

^gx — l^gx Wf^gx 

'Tgx 

Pgy = Pgy + '^0gy 

^gy 



ki 

^Vx = -g sin 9 - ^Vx + Wax 

m 

ki 

^Vy = g COS 9 sin (j) - ^Vy + Wo 


(9) 


where, Wax and Way are WON terms included to account 
for the model imperfections in Q. 

Equations 0 , 0 and 0 together describe the process 
dynamics of the estimator. The resulting system can be rep¬ 
resented as a non-linear function of states, control inputs and 
noise terms. 


X = /(x,u,w) 


Measurement Model 

Observations of the EKE are the measurements from X and 
Y accelerometers, which are aligned respectively with bi and 
62 - Measurement equations can be easily derived from 0, 
after including accelerometer noise terms, which are assumed 
to be Gaussian. 


ki b 

CLx = - V- 

m 

ki b 

Qy = - V. 

m 


V 


Wa 


Wn 


( 10 ) 


where Qx and Oy are respectively the measurements from the 
X and Y axis accelerometers on-board the quadrotor. Here we 
assume that accelerometer biases are random constant values 
which can be compensated for, offline. 


EKF Mechanization Equations 

For the mechanization of the Extended Kalman Filter, the 
discrete state transition matrix Ak should be calculated. For 
this we first calculate F, which is the Jacobian matrix of partial 
derivatives of / with respect to x. Then Ak is calculated by 
discretization of the Jacobian matrix. 


F{t) 


df{x,u,w) 


dx 


ikMk 


Discretization is performed with a truncated Taylor series 
approximation and a sample time of Tg, resulting in. 


Ak — J + F[t)Ts 

In deriving the discrete process noise matrix Qk, we assume 
that noise terms in o and 0 are uncorrelated with each other 
as well as with accelerometer noise terms. 


W — ^gy Wj^gx Wi^gy Wax tCay] 

JK(f) diag \^^gx ^gy Pgx ^Pgy ^ax ^ayl 

Q{t) = G{t)W{t)G^it) 

The first four terms of the (t) are the noise variances of 
gyroscope sensors and their biases. These can be found by 
experimentation with actual sensors. Last two terms, which 
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correspond to the uncertainty in (|^ were approximated first 
and then fine tuned for optimum performance of the estimator. 
Also, 


G{t) 


df(x,u,w) 


dw 


ik,<ik 


Discretization of Q{t) results in Qk- 


Qk= f 
Jo 

Measurement matrix H required for the EKF can be directly 
obtained from (|T0|) as. 



Fig. 5. ARDrone Quadrotor used for experiments 


0 0 0 0 —ki/m 0 

0 0 0 0 0 —ki/m 


Assuming uncorrelated errors in accelerometer measure¬ 
ments, measurement noise matrix Rk becomes diagonal, con¬ 
sisting only of the noise variances of the X and Y accelerom¬ 
eters. 


Rk = diag [cr^ aly] 

For initialisation, all states of the filter are set to zero 
and their error covariances are set to small positive values 
reflecting the uncertainty in initial estimate. With multiple 
experimental runs, it was found that changes of up to 100% 
in the initial values and the noise variances have negligible 
effect on filter performance. We attribute this robustness of the 
estimator to the linear measurement model and not-so-strong 
non-linearities in the process equations. 

EKF state prediction was carried out with the use of a 2"'^ 
order Runga-Kutta integrator. Covariance projection, Kalman 
gain calculation, state update and covariance update equations 
of the estimator take their standard forms as detailed in M- 

ARDrone Quadrotor and the experiments 

The quadrotor platform used for the experiments presented 
in this article is the Parrot ARDrone ifTTll (see Fig. ARDrone 
weighs about 420g including the protective hull and has a 
flight time of about 10 minutes. Straight out of the box, 
ARDrone is an extremely stable quadrotor platform and there¬ 
fore is an excellent platform for quadrotor based research. It 
is equipped with a wide array of sensors including triad of 
accelerometers, triad of gyroscopes, two cameras -one facing 
front and other facing down- and downward pointing sonar 
sensors. All sensor data from the ARDrone are wirelessly 
transmitted to a ground station PC either running Windows 
or Finux. An open source C API is provided which can be 
easily extended to develop application on the ground station 
to process incoming sensor data and to send out control 
commands to the ARDrone. It is also equipped with a pre¬ 
programmed closed source attitude control system, which 
takes care of the low-level stabilisation and control tasks, while 
providing users the ability to develop applications for higher 
level navigational tasks. 

It is desirable to have “ground truth” states trajectories for 
performance evaluation of the proposed estimator. Therefore, 


all our ARDrone experiments were performed in a Vicon 
motion capture environment. The Vicon motion capture system 
uses a set of reflective markers rigidly attached to the quadrotor 
body, which are observed by 8 fixed IR cameras to directly 
compute the attitude and position of the quadrotor with respect 
to the Vicon coordinate frame. 

In a typical experiment, ARDrone was manually piloted 
within the Vicon environment (approximately 6 x 4 x 3 m) 
using a joy stick attached to the ground station computer. 
The inertial sensor data were continuously streamed to the 
ground station computer at 200Hz and were stored for post 
processing. Vicon generated state estimated were also stored 
in a separate PC. Matlab computing environment was used for 
post processing of both inertial and Vicon data. 

A critical parameter that needs to be precomputed for the 
estimator is the rotor drag coefficient Ai. Since a theoretical 
calculation of this parameter is a complex task, we resorted to 
an experimental estimation method. The basic methodology 
adopted here is to obtain the accelerometer measurements 
and ground truth velocity data of a few flight tests. A rough 
estimate of the parameter ki (which incorporates Ai) can then 
be obtained by formulating (lOi as a least-squares problem. 
For the ARDrone, the best estimate for the parameter ki was 
found to be 0.57. This parameter estimation task was run only 
once and the derived ki values was used for all subsequent 
estimation tasks. 


Experimental Results 

During one experiment, the AR Drone was manually oper¬ 
ated within the Vicon environment, moving freely while keep¬ 
ing the height approximately constant. A three-dimensional 
trace of the path taken by the MAV in a typical experiment is 
shown in Fig. The results presented in the following sections 
are based on the data gathered from this experiment. 

Fig. 1^ shows the attitude estimates of the proposed EKF 
together with the ground truth obtained from the Vicon system. 
For comparison purposes, we have also plotted the attitude 
estimates from a generic estimator as detailed in im in 
Fig. [9] It is important to note the improvement in the pitch 
estimate of the proposed estimator over the generic estimator. 
This improvement is more pronounced in places where the 
quadrotor changes its flight direction (for example around 4.6 
and 7.8 sec). During those intervals, the quadrotor undergoes 
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Fig. 6. Three-dimensional flight path of the ARDrone experiment 


high inertial accelerations and the assumption that the ac¬ 
celerometer measurements are dominated by gravitational ac¬ 
celeration fails to hold. Thus generic attitude estimators based 
on this assumption produce erroneous results. As expected 
the proposed EKF attitude estimates agrees more with the 
ground truth because such an assumption is not utilized in 
that design. However, when the quadrotor is not undergoing 
considerable accelerations, the two attitude estimates converge 
and the generic estimator can perform just as well as the 
proposed method. 

Fig. |1 l(a)| and |1 l(b)| present a comparison between the 
errors in the roll and pitch attitude estimates of both the 
proposed FKF and the generic estimator. Fven with the 
proposed FKF, unmodelled dynamics (such as displacement of 
accelerometer from the centre of mass of the quadrotor) causes 
an increase in estimation error when the quadrotor undergoes 
large accelerations. But overall, it is clear that the errors in 
the proposed design are considerably less than those of the 
generic design. 

Fig. [^presents the velocity estimate from the proposed FKF 
together with the ground truth. Again for comparison. Fig. 


10 shows the velocity estimates in a generic design where. 


velocity is estimated by integrating inertial accelerations cal¬ 
culated by compensating the accelerometer measurements for 
gravity. A comparison between the errors in velocity estimate 
obtained from the proposed estimator and the generic estimator 
is shown in Fig. 11(c) where total velocity error is the sum of 
root square errors of both X and Y axes. What is important to 
note is that the proposed strategy produces velocity estimates 
in which errors do not grow with time, while estimating veloc¬ 
ity through direct integration of accelerations as implemented 
in the conventional design leads to a significant drift. As zero 
velocity updates, that can be used to correct this behaviour 
in land vehicles, are no longer viable with an MAV without 
some deliberate control strategies, this points to a significant 
advantage of the estimator proposed in this article. 


Conclusion 

In this article, we presented a novel state estimator for 
quadrotor MAVs, where clear improvements in estimates stem¬ 
ming from the incorporation of quadrotor specific dynamical 
constraints were demonstrated. Our design is based on an 




Fig. 7. Comparison of ground truth and inertial attitude estimates of AR 
Drone. [(^ Roll angle (())), |(b)| Pitch angle (6) 



(a) 



Fig. 8. Comparison of gr ound truth and inertial velocity estimates of AR 
Drone. [(^ X Velocity (V 3 ;),|(b)|Y Velocity (Vy) 


EKF and is capable of estimating both roll and pitch angles 
of the attitude in addition to X and Y components of the 
body frame translational velocities within a bounded error. 
This estimator is applied to inertial data gathered from real 
world flight experiments. The resulting attitude and velocity 
estimates obtained match closely with the ground truth and 
are drift free. 

Before concluding the discussion on the estimator perfor¬ 
mance, we note that our design by itself is not a perfect 
solution to the problem of quadrotor state estimation. We 
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Fig. 9. Comparison of ground truth and inertial attitude estimates of AR 
Drone, obtained from the generic estimator. |(a) | Roll angle (0), |(b)| Pitch angle 

(ff) 





(c) 


Fig. 11. Estimation errors of both estimator designs. |(a)| Roll angle {cj>) 
estimation error, |(b)| Pitch angle (0) estimation error, |(c)| Total velocity 
estimation error 



(b) 


Fig. 10. Comparison of ground truth an d ine rtial velocity es tima tes of AR 
Drone, obtained from the generic estimator. |(a)| X Velocity (Vi),[(^Y Velocity 
(Vy) 


believe that two key improvements need to be made to our 
design. First, an online estimation of the parameter Ai and 
accelerometer biases will improve estimation accuracy and 
ease the filter design process. Secondly, the estimation ^ angle 
and velocity will improve the autonomy of the quadrotor. 
Our cuiTent research focuses on these improvements. 

In addition, we also expect to fuse the inertial information 


with exteroceptive sensors such as cameras and GPS. The two 
cameras in the ARDrone makes it an ideal platform for visual 
Simultaneous Localisation And Mapping (SLAM). One key 
drawback in employing monocular SLAM for MAVs is the 
unavailability of odometry for scale recovery. Another more 
obscure problem is the alignment of camera with the MAV 
body frame. From a control theoretic perspective, orientation 
of the body frame is what matters and misalignment of camera 
and body frames can lead to poor control performance in 
a SLAM only MAV state estimator. Both these problems 
can be solved by tightly integrating the estimation algorithm 
presented here with a monocular SLAM algorithm. We believe 
this to be an exciting research avenue. 
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